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ABSTRACT 

The problem of a particle trapped in a time dependent meridional magnetic 
mirror field and at the same time subjected to a time dependent perpendicular 
electric field (E • B = 0) is considered. Bounce averaged guiding center theory 
is used to derive expressions for the drift velocity components. These drift ex- 
pressions depend in a straightforward fashion on the structures of the magnetic 
and electric fields and on the particle charge, kinetic energy, longitudinal in- 
variant, and bounce period. In the special case of a static magnetic field, no 
electric field, and J = 0 particles, the drift equations are integrated and the 
bounce averaged guiding center trajectory obtained. 
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GUIDING CENTER DRIFTS IN TIME-DEPENDENT 


MERIDIONAL MAGNETIC FIELDS 


INTRODUCTION 

We are concerned in this paper with the guiding center drifts of particles 
trapped in mirror magnetic fields whose field lines are totally contained in con- 
stant longitude— i.e., constant 4 >— planes (in a spherical coordinate system). We 
call magnetic fields of this type meridional magnetic fields, the field of a mag- 
netic dipole being an example. However, no requirement of axi-symmetry is 
imposed on the field. Time variations which preserve the meridional character 
of the field are also allowed. 

Electric fields accompany such magnetic time variations. In addition there 
may exist electrostatic electric fields (VxE = 0). (A source of such electro- 
static fields might be alow-/3 plasma confined by the magnetic geometry.) We 
here assume that the total electric field lies perpendicular to the magnetic field 
(E • B = 0). 

For the field of a static magnetic dipole and in the absence of electric fields, 
Northrop [1966], using guiding center theory, has derived a rather simple ex- 
pression for (4>}y the longitudinal drift rate of a particle averaged over its 
bounce motion (Bounce averages are denoted by ( ).), 



Here c is the speed of light; p. is the magnetic moment of the dipole; r Q is the 
equatorial crossing radius of the guiding center drift shell; and e, W, J, and T 
are respectively the charge, kinetic energy, longitudinal invariant, and bounce 
period of the drifting particle. 


J - m j> V|| ds , (2) 

f ds 

T = K' (3) 

m being the particle mass and the s-integrals being extended over a complete 
guiding center bounce period. Equation (1) for the dipole field drift rate can be 
shown (with some amount of algebra) to be exactly equivalent to the expression 
derived by Hamlin, Karplus, Vik, and Watson [l96l], Northrop's result is 
aesthetically pleasing in that only the familiar particle parameters e, W, J, 
and T appear in it; our results for the drifts in more general fields share this 
property . 

THEORY 

We have found that for other meridional (but not necessarily either axi- 
symmetric or time independent) magnetic fields, the drift rate is simply expressed 
in terms of the particle parameters e, W, J, and T; the quantities a and ^identifying 
the field line upon which the guiding center is instantaneously bouncing; and 
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derivatives of the electric potential <|>. A derivation of these results is com- 
municated here. 

A guiding center drifts across magnetic field lines at a rate which depends 
on its bounce phase. We here assume that the size of the electric field and the 
space and time variations of both the magnetic and electric fields are sufficiently 
small that over a complete bounce period a guiding center drifts only a small 
amount. The bounce path is then nearly periodic, and it is significant to consider 
the rate of guiding center drift averaged over the bounce motion. We are here 
concerned with this average drift of guiding centers. 

The bounce averaged drift rate is resolved into components (a) and (/3) , 
the Euler potentials a and /? being defined in the usual fashion 

B(r, t) = Va(r, t) x V/3(r, t) . (4) 

The drifts (a} and \/3) determine a trajectory in a cartesian two-dimensional 
a, P space, each point of which represents a magnetic field line. The relations 
a = a(r, t), /3 - /3(r, t)may subsequently be used to translate the a, (3 tra- 
jectory into one — from field line to field line — in real space. 

Northrop [l963] provides us with the dynamical equations for averaged 


guiding centers. 



correct to leading order in the adiabatic (m/e) expansion. The quantity K is 
given for each particle by the expression 


mv„ 


K 


+ MB + e 


[«> + tlf ( r ’ t)] 


( 6 ) 


with M the particle magnetic moment (M - Wj/B) . In the absence of time varia- 
tions K is thus the total energy of a particle moving in the magnetic field B through 
the electrostatic potential 4>. 

For meridional magnetic fields we can and shall choose /3 as longitude 4>. 

With this choice V/3 = (r sin d) -1 is totally azimuthal (hence perpendicular 
to the magnetic field) and dfi/d t (r, t) =0. The Euler potential a(r, d, 4>, t) is 
then from Equation (4) any well behaved solution to 


If" - r 2 sin d B f , (7a) 

dd 

- r sin B^j . (7b) 


On the time scale of many bounces K is not a conserved quantity in this time 
dependent situation. 

We now solve Equation (6) for v |( and substitute the result into the defining 
Equation (2) for J. We do so, noting that the infinitesimal arc length ds is 
(from the field line equation) rB dd/B^, dd being the increment in co-latitude 
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subtended by ds. Substituting the ensuing expression for J into Equations (5), 


we obtain 




rB da(MB + e ^) 
2B t? [K - MB - e$ )] x 


(8a) 


0 



rB 

2B t? 


d± CMgj 


e<D) 


[K - MB " eO] 


1/2 I 


(8b) 


The integrals in Equations (8) are carried out at constant a and between the 
mirror co-latitudes, and of the bouncing particle. The multiplicative 
factor of 2 in each appears to compensate for the fact that the d-integral repre- 
sents only half the bounce motion. No contribution results from differentiation 
of the limits of integration since v (1 vanishes at both and d ;: . 

We now rearrange the integrands of Equations (8) 


<*>> 



[K - MB - e$] 


1/2 


B d J_ / rB\ 1 dB 
rB da ^B^ J + 2B da 


\ [k -mb - ed>] _1/2 - 


(K~ eO) 


B 


dB 

da 


+ e 



(9a) 
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The general expression for the electric field [Northrop, 1963] 


E 


1 

c 



( 10 ) 


reduces for the choice fi - <f> to 


E 


e - - <» 17 - V* 

cr sm <9t ^ 


(ii) 


In our model E • B = 0 and it follows that magnetic field lines are equipotentials, 
so that <9$ /da, and d$/d4> remain constant for the integrations occurring in 
Equations (9). Furthermore, to the order in m/e which we are working the 
variation in K over the bounce motion is negligible [Northrop, 1963]. 

Consider now mirror-type, meridional magnetic fields having the properties 

d (rB\ 1 dB 

rB da \B#) + 2B da ~ a o O* (12a) 


1 dB 

B da ~ a i O’ <*>' 1) ’ (12b) 
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a Q and a x being functions independent of distance along the field line a, 0. Fields 
with the properties Equations (12) [and similarly fields satisfying Equations (14) 
following] are selected for consideration because of the simplicity of the resultant 
expressions for the drifts. For fields of the form Equation (12), ( 0 ) becomes 


<0> 


c 

eT 


Ja, 


-t(' 


Wa , + e 



(13) 


In obtaining Equation (13) from Equation (9a) the definitions of J and T were used 


and the particle kinetic energy W - K- e4> was identified. 


Similarly, for magnetic fields with the properties 


B d d /rB\ 1 dB 

rB a0 \B,j] + 2B d4> ~ b o ( a - 1 ) - ( 14a ) 


i aB _ , 

b aa> °i ( a ’ c p < t ) ’ 


(14b) 


the drift ( assumes the form 


(i> 


c 

eT 


J b o-T 


(wb t + e 


a<iA 
<90/ ' 


(15) 


When both the meridional magnetic field and the potential are axi-symmetric, 

(a) = 0. 
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Mirror -type meridional magnetic fields of the form 


B r " ^ (P^“2) g sfn $ w [ sin d f (^- t )] - ( 16a ) 

= -jj f(i>, t ) g(0, t) , (16b) 

= 0 , (16c) 

separable in the coordinates r, d, and 4>, have the properties indicated in Equa- 
tions (12) and (14). In Equations (16) p is an arbitrary number greater than 2, 
and f and g are subject to the restriction that the magnetic field be well defined 
for all r > 0. The axi-symmetric multipole fields are special, curl -free cases 
of Equations (16) and correspond to integral p(>3) , g independent of 4>, and 
f(d) - - dP p _ 2 (cos d )/dd, Pp_ 2 being the Legendre polynomial of order p- 2. 

Equations (16) describe a mirror geometry in a region of space provided that 
at least one magnetic minimum occurs on field lines, i.e., <9B/dd| a ^ t = 0» 
d 2 B/<9d 2 1 a ^ t > 0 for at least one point on B-lines in that spatial region. 

For the field, Equation (16), we choose 

r 2 - p 

a ~ p _ 2 sin d f(d> t ) g(4>, t ) . (17) 

With this choice all gauge freedom in this problem is eliminated, for d<P/da and 

must now be such as to yield the correct electric field through Equation (11). 
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It may now be directly verified that 


a 


o 


2a ’ 


(18a) 


= P 1 

a i p - 2 a 


(18b) 


b Q = 0, 


(18c) 


and 


- 2 1 dg 

l 2 -p g d4> 


(18d) 


Substituting into Equations (13) and (15), we find for the drift components 



(19) 


£ ( 2W 1 <?g 
e \2 - p g d<t> + e ) 


( 20 ) 


Northrop’s dipole field result, Equation (1), is recovered from Equation (19) by 
the identification $ = 0, p = 3, g = -p, f = sind, and (from Equation 17) 
a = - jj. sin 2 d/r = * fJ./ r Q . 

In the general, time dependent case, the complicated a, <fi, t dependence of 
the right sides of Equations (19) and (20) prohibits analytic determination of 
the guiding center trajectory. In simplification we therefore focus attention on 
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the situation where the magnetic field is static and there are no electric fields. 
Now the particle kinetic energy W is constant, <I> is zero, and we obtain 


da 

dd> 


(») 

<*> 


2WT a dg Ij _ pWTl ~ 1 

2 - p g d 2 p~2 


( 21 ) 


by dividing Equation (20) by Equation (19). 

Even in the form Equation (21), the a,<£ dependence of the bounce period T 
inhibits further progress. The case of J = 0 particles is, however, one limiting 
situation which can be handled. These particles have no component of motion 
along the magnetic field and drift so as to always remain at a local magnetic 
minimum, dB/dd\ a 0. The bounce period T for such particles defined in the 
J-*0 limit is non-vanishing. For J = 0 particles, Equation (21) reduces to 


da _ 2a dg 
dcp ~ pg dcfi 


( 22 ) 


Equation (22) can be integrated to yield 

a = Ag 2/p , (23) 

A being a constant determined for each trajectory . In order that the magnetic 
field be well defined for all r > 0, g must be a periodic function of <t>. It follows 
then from Equation (23) that a particle executes closed loops in a, /3 space, 
indicative of the fact that after each real space excursion of 2^ in longitude a 
particle returns to the same field line for this static situation. 
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The time development of the trajectory may now be traced using Equation 
(23) coupled with Equation (19) simplified for static B, 4> = 0, and J = 0, 

d 4 > c p c p , 

HC = 2/t - (24) 

Integration of Equation (24) depends, of course, on a choice of g, and for most g's the 
4 > -integral cannot be explicitly evaluated. For the interesting case of small azimuthal 
modulations of the form g = 1 + e cos" 4 >, e « 1, n > o, it is possible, however, 
to develop from Equation (24) the time history of 4 > as an asymptotic series in e . 
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